Many-body Floquet theory of laser-induced phase transition in quantum magnets 
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We theoretically predict a non-equilibrium phase transition in quantum spin systems induced by 
laser, which provides a purely quantum-mechanical way for a coherent control of magnetization. 
Namely, when a circularly polarized laser is applied to a spin system, the magnetic component of 
laser induces a magnetization normal to the plane of polarization, leading to an ultrafast phase 
transition. We first demonstrate this numerically for an S — 1 antiferromagnetic Heisenberg spin 
chain, where the topologically-ordered Haldane state gives way to a new state with an emerging 
magnetization. We then elucidate its physical mechanism by applying a many-body Floquet picture. 
The theory also indicates that the phenomenon should occur in general quantum spin systems with 
a magnetic anisotropy. The required laser frequency is in the THz range, and the intensity within 
a prospective experimental feasibility. 

PACS numbers: 42.55.Ah, 75.10.Jm, 75.10.Pq, 75.30.Gw 
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There is a growing fascination with physics of non- 
equilibrium systems, which is becoming an important 
topic in condensed-matter and other fields of physics. For 
electron systems, a host of novel phenomena induced in 
non-equilibrium situations, such as photo-induced Mott 
transitions Q , photo-induced topological transitions [|| , 
etc, have been fathomed. It is a hot topic in cold 
atoms as well, where quantum simulation is being re- 
alized ||. Now we pose a question: "can we propose 
novel non-equilibrium phenomena for quantum spin sys- 
tems?" Spin systems, a basis for the physics of magnets, 
are many-body in nature and have various applications 
in engineering and, more recently, in spintronics. This 
paper proposes a novel, non-equilibrium way to coher- 
ently control many-body states in spin systems with laser. 
This is in contrast to the control of single quantum states, 
e.g., qubits, which is becoming important in the field of 
quantum computation Q . Instead, we want to develop a 
method to control the collective phenomena, e.g., phase 
transitions, with laser. 

Interaction between matter and laser, the most pow- 
erful tool for coherent control, has been mainly used in 
electron and atomic systems. So the first question one 
might ask is " can spin systems be manipulated by laser 
as well?" The second question, no less important, is "can 
it be done coherently on a many-body state?" Control of 
spin systems by lasers is a hot topic due to recent exper- 
imental advances [|[ The key of our study is to use 
the magnetic component of circularly polarized lasers with 
photon energy far below the electron energy scale. It was 
demonstrated recently that magnetic fields from lasers 
in the THz regime can directly access to spin dynamics 
without disturbing the charge degrees of freedom || . The 
merit is that incoherent processes resulting from charge 
excitations can be prevented. While the setup of Ref. || 
was within the linear response regime, here we extend the 
avenue to the non-perturbative regime by proposing that 



one can realize a "laser-induced phase transition" with 
perfect quantum coherence. We show that, with circu- 
larly polarized laser, we can induce net magnetization 
in a quantum antifcrromagnet. First, it is demonstrated 
numerically with the infinite time-evolving block decima- 
tion (iTEBD) [[jj |j| for one-dimensional spin chains. The 
motion of the induced magnetization can be controlled 
as desired by shaping the laser pulse, hence provides a 
new technique for coherent spin manipulation. 

The dynamical phase transition requires a theoretical 
treatment which goes beyond the linear response theory. 
Floquet theory is becoming a standard theoretical pic- 
ture for studying quantum systems under periodic driv- 
ing ^ ^-Ejj. Namely, the time-periodicity enables us 
to cast a time dependent problem into a static effective 
model governed by the Floquet Hamiltonian. This has 
proved to be a powerful method in the theory of "Floquet 
topological insulators" B pL |o|, and has also been ap- 
plied to many-body problems such as the photo-induced 
We usually have to deal with a 
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Mott transition 
Hamiltonian in an infinite-dimensional Floquet matrix 
form as a penalty for going from the time-dependent 
problem to a time-independent one. We find here, how- 
ever, a dramatic reduction of the dimension for the Flo- 
quet Hamiltonian due to a symmetry by applying the 
many-body Floquet method to spin systems in circularly 
polarized laser. The reduced model, which we call the "ir- 
reducible Floquet model" , turns out to be a spin system 
with a static slanted magnetic field. The emergence of 
the laser-induced magnetization can be understood with 
this static model, which we confirm by comparing exact 
diagonalization results with iTEBD. Since the above dis- 
cussion is applicable also for more than one dimension, 
the induction of magnetization is independent of the di- 
mension. 

While the magnetization can be induced in any di- 
mensions, one-dimensional antiferromagnet has, quan- 
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turn mechanically, a special interest. The system is 
known to be gapless if the size of each spin, S, is a half- 
odd-integer, and gapped if S is an integer. This is the cel- 
ebrated Haldane's conjecture 12| now established by 
intensive analytic, numerical, and experimental studies. 
Specifically, the ground state (GS) of S = 1 spin chain, 
known as the Haldane phase, is topologically protected 
by a symmetry [ H, a nd is characterized by a "string 
order parameter." JL5| This enables us to demonstrate 
that the laser-induced non-equilibrium phenomenon we 
propose here is indeed a phase transition from a state 
having a topological order into a state with finite magne- 
tization which replaces the topological order. As a fur- 
ther indication that the present transition is genuinely 
a non-equilibrium process, we show that the necessary 
strength of the magnetic field can be smaller than the 
static excitation (Haldane) gap, which is not the case 
with conventional magnetization processes in static mag- 
netic fields. 

Let us begin by introducing the model used in the 
present study. We consider a one-dimensional Heisen- 
berg antiferromagnet with a Hamiltonian, 

u = jj2 s i- s j+ D J2( s if ( J> °)' w 
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FIG. 1: (color online), (a) A quantum spin system applied by 
a circularly polarized laser is schematically depicted, (b) A 
typical evolution of the magnetization induced by the laser, 
as calculated with iTEBD for D/J = 0.25, A/ J = 1.0, and 
Q./ J — 1.4. The range of the plot corresponds to about two 
cycles of the laser (a time interval < t < 9. 6 J -1 ), (c) 
Time-evolution of magnetization, M z , along z-axis, and phase 
difference, AO, in the xy-p\a,ne between the magnetic field of 
laser and magnetization. 



where the term proportional to D is a magnetic 
anisotropy of a single-ion type. The model is sim- 
ple, and also realistic for describing, e.g., organic com- 
pounds such as Ni(C 2 H 8 N2)2N0 2 C104 (NENP) and 
Ni(C 5 D 14 N2)2N3(PF 6 ) (NDMAP). D/J in NENP and 
NDMAP is estimated from experiments to be about 
0.18 H and 0.25 ffj, respectively. Here we set D/J = 
0.25. 

We consider a sudden switch on of the laser at t = 
starting from the GS of (Q) which is in the Haldane 
phase jl8|. The state evolves according to the Hamil- 
tonian, 



H(t) = H - A( e - [nt s 1 
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for t > 0, where 2A (= B) and fl are the amplitude 
and frequency (photon energy) of the laser, respectively, 
and S^ ot = J2i the raising and lowering operators 
for the total spin. Here we assume that only the mag- 
netic component of the laser couples to the system, and 
that the laser is applied from the direction of z-axis. 
Thus each spin feels a magnetic field rotating in the 
xy-plane, B(S^ ot cos fit + S% ot sin fit), which is expressed 
above in terms of the spin raising and lowering operators 
(S^t = S* ot ± iStot)- The time evolution of the induced 
magnetization M (= (Si) here) calculated with iTEBD 
is displayed in Fig. |l|. The magnetization starts from 
zero and grows in magnitude with time in a processional 
movement. Remarkably, M tends to point in the z di- 
rection, despite the external magnetic field being entirely 
in the xy-plane. 



There are two conditions for the emergence of mag- 
netization. One is the anisotropy of the spin system. 
In a SU(2) symmetric system, the spin dynamics be- 
comes trivial since Ho commutes with S^~ ot , S^ ot and 
Stot = J2i I n tne present system, the single- ion mag- 
netic anisotropy, i.e., the D term, lifts this constraint. An 
Ising-type anisotropy (XXZ model) has a similar effect as 
well. The second condition is to use a circular polariza- 
tion. In linear polarization, emergence of magnetization 
perpendicular to an applied magnetic field is prohibited 
by spin-inversion symmetry. Circular polarization breaks 
the symmetry and induces magnetization initially point- 
ing in the xy-plane, and then slowly moving toward the 
z-direction. As for the precessional movement, there is an 
intuitive way to grasp the behavior with a semi-classical 
spin picture. Fig. |l](c) plots the evolution of M z , which is 
non-monotonic. If we compare this with the relative an- 
gle between the directions of the external magnetic field 
and magnetization defined by Ad = fit— arctan(M J, /M a: ) 
mod 7T, the increase and decrease of M 2 corresponds to 
the positive and negative AO regions. This indicates that 
M = (M x , M v , M z ) follows a semi-classical equation, 
M = 7M x B(t), where B(t) = B(cos fit, smflt, 0) is 
the magnetic field of the laser with 7 a positive constant. 
If we use fields with opposite polarization (left vs. right), 
the magnetization points in the other direction. 

Although we have an intuitive understanding on why 
the magnetization points in the z-direction, the emer- 
gence of the magnetization itself is purely a quantum 
process. In order to clarify the mechanism, we plot how 
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FIG. 2: (color online), (a) ^-dependence of M z at various 
times t. (b) ^-dependence of M z for t/J' 1 = 2.4. The solid 
line is a fitting to M z oc [A/ J) 2 . 



the induced magnetization depends on the laser ampli- 
tude and frequency in Fig. ||. From the oc (A/ J) 2 depen- 
dence for small A, we conclude that this is a second-order 
nonlinear process in terms of the external magnetic field. 
Figure ||(b) shows the ^-dependence of M z at various 
times t. As time advances, the peak in M z develops 
around Q/ J ~ 1.4, implying a resonance around this fre- 
quency. We note that this energy scale is an order of 
magnitude larger than the Haldane gap A/ J ~ 0.26 J, 
which is a hallmark of the contribution from high energy 
excited states. In the following, we explain the resonance 
behavior with the help of the Floquet theory combined 
with exact diagonalization p9[ . 

The Floquet theory is a temporal analogue of the 
Bloch theorem and maps a system with a time-periodic 
Hamiltonian H(t) to a static problem. This is done 
through discrete Fourier transformation of the Hamil- 
tonian H(t) = ^ m e~ lmQjt H mi and the time-dependent 
Schrodinger equation is casted into an eigenvalue equa- 
tion 



{H n . 



^5 mn )\<S> m ) = 



(3) 



where e is the Floquet quasi-energy level. The solution 
of the original Schodinger equation is given as |^(t)) = 
^2 m e - 1 ( E +« in ) t |$ m ). I n many-body Floquet theory, we 
can think the Fourier index m as a flavor or pseudo- 
spin which runs from — oo to +oo and the term mQS mn 
then acts as an effective "flavor Zeeman term". Now, 
the problem becomes time- independent, but a price to 
pay is that the size of the Hilbert space becomes infinite 
times larger due to this flavor index. Fortunately enough, 
however, in spin systems in circularly polarized light, the 
inter-flavor coupling in the Floquet model is given by 



-ASl 



-AS: 



i.e., the raising (lowering) operator of the flavor index 
is proportional to the raising (lowering) operator of the 
total spin, and [H ,S^ ot ] = is satisfied. Therefore, the 
quantity (S* ot — m) is a good quantum number. Thus the 
Floquet Hamiltonian Eq. (|6|) becomes block-diagonal in 



terms of this quantum number and each block is governed 
by what we call the "irreducible Floquet Hamiltonian" , 



H It = J^Si-Sj+D^iSf)* 
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where B = 2A, and the z-direction magnetic field f2 
comes from the "flavor Zccman term". Different blocks 
are related through a simple energy shift in an integer 
multiple of tt. The obtained effective Hamiltonian, rep- 
resents nothing but a spin system in a slanted magnetic 
field B e ft = (B, 0, fl) (Fig. ||(a)). The above discussion is 
applicable to any lattice structures other than the chains 
considered here. We can note in passing that the effective 
description in the Floquet picture has a resemblance with 
the electron spin resonance (e.g., Ref. 20), where the role 
of the external magnetic field is now played by the laser 
frequency f2. Now it is not difficult to understand the 
mechanism of the laser-induced magnetization process. 
Excited states with S£ ot — ±j (j > 0) are initially de- 
generate due to the spin inversion symmetry. When the 
laser is applied, the "longitudinal magnetic field" ft lifts 
the degeneracy due to the Zeeman effect and the Floquet 
energy levels split as Ej — !• Ej TJ^ as shown in Fig. ||(b). 
The "transverse field" B acts to hybridize the states hav- 
ing different Scot's- The hybridization becomes larger 
as the Floquet quasi-energy levels come closer. Since 



+j states approaches the GS while 



-J 



states departs from the GS, the GS primarily hybridizes 
with S* ot = +j states, which is precisely why a positive 
net magnetization appears. 

It is non-trivial to take the thermodynamic limit since 
we have to consider the hybridization with infinitely 
many levels. To clarify this, we introduce the "magnetic 
density of states" (MDOS), defined as MfS(Ei - E), 
which enables us to capture the behavior of the Zeeman 
splitting for the entire many-body states. MDOS for var- 
ious values of Q calculated with diagonalization of the 
irreducible Floquet Hamiltonian is shown in Fig. ||(c) for 
N = 8 and A/ J = 0.1, where it is normalized by 2^. 
MDOS, which is zero at fi = due to the S% ot = ±j 
degeneracy, is seen to develop positive (negative) peaks, 
which move to low (high) energies with ft. When the 
positive peak overlaps with the GS, the hybridization be- 
comes most prominent, and thus, leads to the resonance- 
like behavior. The Floquet theory is so powerful that one 
can reproduce the iTEBD result of the time evolution of 
M z . Figure ||(d) plots M z against Q at various times 
calculated by combining the Floquet theory and the ex- 
act diagonalization for N — 8 E^). We can see that the 
result agrees well with the previous one (Fig. f§(b)). The 
peak is located around ft/ J ~ 1.9, and the deviation 
from the iTEBD result (peaked around Q./J ~ 1.4) is 
attributed to a finite-size effect. 

Finally, we show that the current phenomena is ac- 
tually a laser-induced phase transition from Haldane's 
phase to the magnetized phase. To demonstrate this ex- 
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FIG. 3: (color online), (a) The Floquet theory maps the spin 
system in circularly polarized laser with amplitude B and 
frequency SI to a static spin model with a slanted magnetic 
field B c ff — (B, 0, SI), (b) An effective description of the laser- 
induced magnetization, where the z-direction magnetic field 
SI shifts the many-body energy levels by "Zeeman splitting" of 
the Floquet energy spectrum, (c) Magnetic density of states 
(MDOS) in units of states/J per site for SI/ J — 0.5, 1.5, 
and 2.5. (d) fi-dependence of M z at various times. Both (c) 
and (d) are obtained by exact diagonalization of the N = 8 
irreducible Floquet Hamiltonian Eq. (pi). 



plicitly, we calculate the string order correlation |15j with 
iTEBD. Namely, Haldane's phase, being a topologically 
protected phase, has no local order parameter, but is in- 
stead characterized by the string order parameter, which 
is defined as lim I ._ >00 O str (r) from the string-order corre- 

lation function, O str (r) = (S§ exp (itt^^ 1 S z ^jS z ) fjj|. 
The string order parameter can capture the Haldane 
phase, which can be considered as a superposition of 
states in which an arbitrary number of \Sf = 0)'s are 
inserted between |1) and | — 1) in a Neel order in a Sf- 
basis, e.g., | . . . , 1, 0, 0, 0, —1, 0, 1, 0, 0, — 1, . . .). In the ini- 
tial GS, we are in the Haldane phase with the string or- 
der parameter indeed finite as seen in the t = result of 
Fig. ||(a). For t > 0, the string order correlation O st i(r) 
begins to decay exponentially with distance r. This in- 
dicates that the laser-induced phase transition from the 
Haldane phase to a magnetized phase happens as soon 
as we have a finite magnetization. Indeed, we find that 
the decay of the string order has a direct (and even an- 
alytic) relation with the uprise of the z-component of 
magnetization. We can in fact quantitatively correlate 
the magnetization with the string correlation length £ 
(defined by fitting O stI (r) cx exp(— r/£)). Figure |](b) 
plots £ as a function of M z . We can deduce an analytic 
form for the relation between £ and M z as follows. 

A flip of a single spin from |0) to |1) or from 
| — 1) to |0) in the Haldane phase gives a fac- 
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FIG. 4: (color online), (a) The string order correlation func- 
tion calculated with iTEBD at various times, (b) The inverse 
of string correlation length, l/£, plotted as a function of M z . 
The result of iTEBD (circles) and the analytical prediction, 
= -ln(l - 2M Z ) (line) are shown. 



tor —1 in the string order correlation function, 
and acts like a "disorder" in the string-ordered 
states (e.g., | . . . , 1, 0,0, 0, -1,0, 1,0, 0, -1, . . .) ->• 
...,1,0, "1", 0,-1, 0,1, 0,0,-1,...)). The probability 
of having such a disorder is M z per site, and by estimat- 
ing the probability of having k disordered sites out of r 
sites, we are led to an expression, 

r 

O st r(r) ~ So £(-l) fc r C k (M z ) k (l - M z ) r - k 
= So(l-2M z ) r , 

where Sq is the initial (i.e., at M z = 0) value of the string 
order. Namely, we have = — ln(l — 2M Z ), with which 
the result in Fig. ||(b) has an excellent agreement in the 
small M z region. 

Finally let us discuss the experimental feasibility. A 
pump-probe experiment should be most realistic, where 
the pump laser is required to be strong and circularly po- 
larized, and the induced magnetization can be measured 
by Kerr rotation. The necessary pump strength depends 
on the sensibility of the probe. We make an estimation 
for NDMAP, whose exchange interaction is J ~ 2.8 meV 
with D/J ~ 0.25 0. While we have developed Floquet 
picture for a continuous laser application, a pulse con- 
taining a few cycles of laser can induce a magnetization 
as illustrated in Fig. [|, and in order to induce M z ~ 0.01 
with few cycles we need A/ J ~ 0.1. This corresponds 
to B = 2A = 0.2J ~ 5.6 T. The optimum frequency 
is, from the resonance frequency, 1.4 J ~ 3.9 meV^ 1 
THz. Currently available intensity of the magnetic field 
in laser is ~ 0.13 T (5| and is still one order smaller than 
the required strength. However, THz laser techniques are 
rapidly advancing and we expect our theoretical proposal 
become experimentally feasible in the near future. 

Finally, let us comment on the outlook of laser-induced 
magnetization. While we have presented the result for 
one-dimensional spin systems, our discussions with the 
Floquet theory do not depend on the dimension nor the 
size of the spin. There are many interesting spin related 
phenomena in higher dimensions, such as the effect of 



frustration and emergence of spin liquid phase, to name 
a few. Moreover, quantum spin systems have interdisci- 
plinary spin-offs to e.g., atomic systems in a cavity. It 
is an interesting future problem to study laser-induced 
phase transitions in such systems, and the theory pre- 
sented here is expected to play an important role. 

We acknowledge illuminating discussions with S. Fu- 
ruya, N. Tsuji, and P. Werner. The computation was 
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Center, Institute for Solid State Physics, the University 
of Tokyo. This work is supported in part by Grants-in- 
Aid from JSPS, Grant No. 23740260. 
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Supplemental Material 

FLOQUET THEORY 

Here we elucidate how the application of Floquet the- 
ory, which usually leads to an infinite-dimensional ma- 
trix equation, to the present quantum spin system can 
be cast into a finite-dimensional form. Floquet theory is 
a mathematical technique to treat time-periodic differen- 
tial equations, which is a temporal analogue of the Bloch 
theorem for spatially periodic systems. When applied to 
the time-dependent Schrodinger equation, 

i^|*(*)>=W(*)|*(t)>, (5) 

with time periodicity H(t + T) = H(i) (T: period), the 
Floquet theorem dictates that the solution should have 
a form \^{t)) = c~ le *|$(£)), a phase factor involving e 
called Floquet quasi-cncrgy times a time-periodic wave 
function (Floquet state) with |$(i + T)) = |$(t)). With 
both %{t) and |$(£)) periodic in t, we can make a discrete 
Fourier transformation, 

m 

!$(*)) = 5~V im °*|$ m ). 



When these are plugged into Eq. (||), the Schrodinger 
equation is casted into a time- independent eigenequation 
in a matrix form, 



J2(H n -m - mn6 mn )\$ m ) = (6) 



which can be thought of as an equation for "photon- 
dressed" (Floquet) modes. 

For the present case of one-dimensional Heisenberg an- 
tiferromagnet with a single-ion anisotropy in a circularly 
polarized field with the Hamiltonian 



H(t) = J^2S i -S 3 +D^2(SD 2 ~A(e- int St ot +e int S^ t ), 



the Oth component (iJ„_ m with n = m) is Hq = 
JY^i j &i ' Sj + DJ2i(Sf) 2 , while H±x just corresponds 
to H±i = —ASf ot . The eigenvalue equation @ then 
simplifies into a tridiagonal form, 











Hn — H-l. 








H 






H+i 

H Q + n 

H-x 






H+x 



■J 



( • \ 

l$ 2 ) 

V i J 



( • N 

|$ 2 ) 

|$°) 

v : y 



(7) 



where f2 represents 57 times an identity matrix. The Flo- 
quet formalism is schematically depicted in Fig. ^|(a), 
where replicas of the original system are prepared for 
different photon-dressed (Floquet) modes m, for which 
the term in the Hamiltonian with a phase factor e T ' m 
induces a transition between the Floquet mode m and 
mil. 

The matrix representation of Eq. (^) is in general 
infinite-dimensional, since m takes from — oo to +oo. In 
the present case, however, S? ot — m is a good quantum 
number because the term S^ t appears in the Hamilto- 
nian with a phase factor e T ' ot , so that the term simul- 
taneously changes the Floquet index and 5 t z ot by ±1, re- 



spectively, as shown in Fig. ||(a). This implies that the 
Floquet matrix can be put into a block-diagonal form as 
shown in Fig. [|(b). We call the Hamiltonian that acts 
within the blocks an "irreducible Floquet Hamiltonian." 
If the system size (i.e., the total number of spins) is finite, 
each block is finite-dimensional since S* ot is bounded as 
— N < S* ot < N. Thus, even in the presence of the 
time-dependent external field, we can readily solve the 
eigenvalue equation with an exact diagonalization as far 
as finite systems are concerned. 

Interestingly, the irreducible Floquet Hamiltonian is 
formally equivalent to the Hamiltonian for an S — 1 chain 
with longitudinal and transverse magnetic fields. Since 
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FIG. 5: (a) Structure of the Floquet Hamiltonian for a spin system in a circularly-polarized light field, (b) Block-diagonal 
structure of the Floquet Hamiltonian reshuffled in terms of the good quantum number 

'-'tot 



the 5* t o t term connects the sectors which have 5t ot differ- 
ing by ±1, B(= 2A) acts as the transverse magnetic field. 
On the other hand, f2 act as the longitudinal magnetic 
field. We can see this because the matrix components 
in the same S* ot sector are Hq — mO, which translates 
to Hq — QSi ot since S z ot — m is constant within each 
irreducible Floquet matrix (where the constant, being ir- 
relevant, can be set to 0). We end up with an effective 
Hamiltonian, 

fflr.Fl. = jY, S ^ S J+ D Y,( S ^ 2 - nS tot~ BS tof (8) 

The time evolution of M z can then be calculated 
from the exact diagonalization result for the eigenval- 
ues ({e Q }) and eigenvectors ({|$q}}) of Eq. (||). We can 
reconstruct the solution of the original time-dependent 



Schrodinger equation as \^f(t)) = Yl a c al*a(*))) where 
I = Em e " i(eQ+m0) *l*™)> and the coefficient c a , 
with \c a \ 2 = 1 for normalization, is determined 
from the initial condition. For example, in the case in 
which the application of laser begins suddenly at t = 0, 
c a = (^ a {t = 0)1*0)1 where |*o) is the initial state, i.e., 
the ground state (GS) of H = JJ2i,j SrSj+D^iiSf) 2 - 
The magnetization per spin M z then evolves as 

+ l4c Q e i ^-^ t (<S>f\SUN\<S>™)+h.c.}, 

where the first (second) term is the ^-dependent (t- 
independent) part. 



